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Abstract. We report an X-ray study of the composite supernova remnant G327.1-1.1, with particular emphasis on 
its thermal emission. By virtue of a combined spatial and spectral analysis, we have been able to model the X-ray 
emission of the remnant as a sum of two components: a non-thermal component, due to the pulsar nebula and 
the pulsar itself, and a thermal component, of which we have analysed spectrum and morphology, after proper 
subtraction of the plerion. We discuss three possible interpretations of the thermal emission of G327.1-1.1: pure 
Sedov expansion, expansion through a inhomogeneous medium with evaporation from ISM clouds, and radiative 
expansion. On the light of our new data and interpretation, we have re-derived all the physical parameters of 
this SNR. In the framework of Sedov or radiative expansion we derive a longer age than previously estimated 
(1.1 x 10 4 ), thus not requiring a high velocity for the pulsar. 

Key words. Acceleration of particles; Radiation mechanisms: non-thermal; ISM: clouds; ISM: individual object: 
G327.1 1.1; ISM: supernova remnants 


1. Introduction 

The class of supernova remnants (SNRs) contains a wide 
variety of objects, that cannot be accounted for by the tra¬ 
ditional subdivision in just two categories, namely shell¬ 
like remnants (powered by a blast wave expanding in the 
ambient medium) and filled-center (or Crab-like, or pie¬ 


showing that a centrally peaked (thermal) X-ray emission 
characterizes SNRs in which the blast wave is expand¬ 
ing in a cloudy interstellar medium and induces a delayed 
evaporation of the clouds. This model has been applied to 
several SNRs, although its validity is still under debate. 
For instancepdho et ah|([1994 ) found that the radial profiles 


rionic) remnan t s (powered by a spinning-down pulsar). 
The class of composite SNRs has m fact been introduced 
(Helfand & Becker 1987] ) in order to arrange a number of 
objects that neither could be classified as pure shell-like 
nor as pure filled-center remnants. This classification was 
originally only morphological, by requiring just the coex¬ 
istence of a shell and of a center-filled component: the un¬ 
derlying belief was that the inner component always has 
a plerionic nature. Such scenario, originally devised for 
objects in which both components are detected in radio, 
seemed to apply also to those appearing as pure shells in 
radio, but with a centrally peaked emission in X rays. A 
natural explanation was that in these SNRs the central X- 
ray component is due to synchrotron emission of plerionic 
nature, which happens to be too weak in radio. 

From the analysis of the X-ray spectra however it be¬ 
came clear that for a large fraction of composite SNRs 
(the so-called “mixed-morphology” remnants) the emis¬ 
sion from the central component is thermal. Along this line 
(White fe Long| 1991, hereafter WL) presented a model 


of SNR W44 are nicely fitted by a WL model; but Harru: 


et al. (1997) then showed that this model fails in dating 
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the SNR, and suggested an alternative scenario based on 
the idea that the shock has reached the radiative phase. A 
more detai led radiative model for W44 has be en recently 
presented ( Cox et al. 1999 : Shelton et al. 1999 ), which in¬ 
cludes thermal conduction and a density gr adient in the 
ambient medium. An alternative suggestion (Petruk 2001) 
is that even a shell-type SNR which expands in a strong 
density gradient may mimic a center-filled morphology, if 
seen with a peculiar orientation. Further investigation is 
still required to clarify which scenario is most appropriate 
to describe the majority of thermal X-ray SNRs with a 
center-filled morphology. 

On the observational side, the diagnostics and even 
the recognition of mixed-morphology SNRs are not easy 
tasks. The determination of the nature (thermal or non- 
thermal) of the central X-ray component on the basis of its 
X-ray spectrum is in principle straightforward. However 
a thermal and a non-thermal component may coexist in 
the same remnant, in which case both a high spatial res¬ 
olution and a broad-spectrum sensitivity are required in 
order to study them separately. The coexistence, in com- 
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posite SNRs, of more components typically implies a dy¬ 
namical interaction between them, a study of which may 
provide important clues on the physical conditions. For 
instance, an X-ray spectral determination of the pressure 
in the shell component, together with the requirement of 
a pressure balance between the thermal remnant and the 
plerion, allows an estimate of the plerionic magnetic field. 

Although there are various composite SNRs that in X 
rays show both a shell-like and a plerionic component, up 
to now only in one case (W44) a plerionic component is 
clearly seen coexisting with a centrally peaked thermal X- 
ray component. In this paper we suggest that also SNR 
G327.1 1.1 may present such coexistence, and be in vari¬ 
ous respects similar to W44. 

G327.1 1.1 is a composite SNR with a rather complex 
structure. Details of its morphology have been revealed 
only recently at radio wavelengths , by th e 843 MHz MOST 
SNR survey ( Whiteoak fe Green 1996 ). In particular its 
plerionic component is substantially offset with respect to 
the centroid of an otherwise rather symmetric remnant. 
A recent Chandra observation has in fact confirmed the 
presence of a compact source (P.O. Slane, private commu¬ 
nication) on the tip of a finger-like feature extending from 
the radio plerion to the north-west direction, which also 
coincides with a hard X-ray source (Sun et al. 1999, here¬ 
after SWC). On the other hand, the thermal component of 
G327.1 1.1 has not been mapped in detail, and its origin 
is still uncertain. SWC interpreted it as a thermal shell, 
and applied the Sedov solution to derive the shell param¬ 
eters. However, the presence of the strong non-thermal 
nebula prevented them from showing the putative limb- 
brightened X-ray morphology of the thermal emission. In 
this work, we attempt to derive the spatial morphology of 
the thermal component. 

The plan of the paper is the following: in Sect. 2 the 
observation is described; Sect. 3 presents results from a 
joined spectral and spatial analysis of the data; Sect. 4 
is devoted to a discussion on the nature of G327.1 1.1; 
Sect. 5 concludes. 


2. Observations 

The BeppoSAX imaging instruments are the Medium- 
Energy Concentrato r Spectrometer (MECS; 1.8-10 keV; 
Boclla et al. |1997) and the Low-Energy Concentrator 


Spectrometer (LECS; 0.1-10 keV; Parmar et al. 1997). 


The MECS consists of two grazing incidence telescopes 
with imaging gas scintillation proportional counters in 
their focal planes. The LECS uses an identical concen¬ 
trator system as the MECS, but utilizes an ultra-thin en¬ 
trance window and a driftless configuration to extend the 
low-energy response to 0.1 keV. The fields of view (FOV) 
of the LECS and MECS are circular with diameters of 37' 
and 56', respectively. In the overlapping energy range, the 
space resolution of both instruments is similar and cor¬ 
responds to 90% encircled energy within a radius of 2(5 
at 1.5 keV. In addition, the BeppoSAX payload includes 
two high energy instruments - the High Pressure Gas 


Scintillation Proportional Counter (HPGSPC; 5-120 keV; 


|Manzo et al. 

1997) 

and the Phoswich Detection System 

(PDS; 15-300 keV; 

Frontera et al. 1997). 


The region of sky containing G327.1-1.1 was observed 
between 1999 March 14 and March 16. The pointing di¬ 
rection was R.A.=15 h 54 m 30 s , Dec.=-55 o 04 , 36" (J2000). 


In order to avoid solar scattered emission and other con¬ 
taminating effects, data were limited to intervals when the 
elevation angle above the Earth’s limb was > 4° and when 
the instrument configurations were nominal. The screened 
exposures in the LECS and MECS are 31.0 ks and 82.5 ks, 
respectively. The LECS and MECS images are shown in 
Fig-0- G327.1 1.1 is located, in both images, near the field 
center, as a resolved peak centered near R.A.=15 h 54 m 29 s , 
Dec.=—55°03'48 ,/ . (J2000). 

An unresolved peak not related to G327.1-1.1, located 
at position R.A.=15 h 56 m 09!6, Dec.=-55°03Tl", with a 
hard power-law spectrum, is also present in MECS2 data. 
However it neither appears in LECS nor in MECS3 data, 
nor it is present in existing catalogs (e.g. RASS). Since 
the spectral response of MECS2 and MECS3 are similar, 
if this source were real the only way to miss it in MECS3 
data would be to have it right behind a rib of the strong¬ 
back: however only soft sources can be hidden completely 
in this way; while ribs are transparent above 5 keV. Of 
course it cannot be an effect of source variability, because 
the integrations of the two MECS are virtually contem¬ 
poraneous. Moreover in a nearby archival BeppoSAX ob¬ 
servation (Target: HD141926, Obs.Code: 20309004), this 
source again appears in MECS2, while it does not in 
MECS3. The conclusion we have reached is that the peak 
seen in MECS2 is a ghost image. It cannot be due to 
problems in the electronics, because it is also visible in 
the archival MECS2 image, while tests on dark and bright 
earth exposure close in time to our observation do not re¬ 
veal anything abnormal. The peak is presumably produced 
by a bright source located off-axis (but that has not been 
identified), and therefore we have simply not considered 
it in the further analysis. Anyway it does not seem to af¬ 
fect the data in the direction of G327.1 1.1. We thank 
F. Fiore and G. Cusumano, of the BeppoSAX team, for 
their valuable advice during the analysis of the nature of 
this source. 

3. Results 

3.1. Spectral analysis 

The thermal and non-thermal components of G327.1 1.1 
are only partially resolved by BeppoSAX, and therefore an 
accurate analysis is required to separate them. This anal¬ 
ysis should take advantage of the differences in the spectra 
of the (soft) thermal and (hard) non-thermal components, 
as well as of the differences in their position and size. 

We have performed spectral fits using different spatial 
and spectral selections of the available photons, and com¬ 
bined the results. In all cases we exclude photons with 
energies which are outside the ranges 0.1-8 keV for LECS 
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LECS 0.1-2.0 keV 



MECS 2.0-10.0 keV 



Fig. 1. 0.1-2 keV LECS image (top) and 2-10 keV MECS 
image (bottom) of the region of G327.1 1.1. Pixel size is 
64" and 32" respectively, and a smoothing of 2 pixels has 
been applied. The images are exposure, vignetting and 
background corrected. LECS contours are at 0.7, 1.3 and 
1.8 x 10 -4 cnt pix -1 , while MECS contours are at 1/16, 
1/8, 1/4 and 1/2 of the peak value (1.2 x 10 -3 cnt pix -1 ). 


and 1.6-10.5 keV for MECS. Fittings have used, sepa¬ 
rately or in combination, the following two emission com¬ 
ponents: a non-thermal component, modelled as a power 
law with photon index 7 , and a thermal one, modelled 
as due to ionized plasma hit by the shock wave (mekal 
model in XSPEC V11.0). Reference abundances are from 
Anders fe Grevesse| 1989, while the interstellar absorp¬ 


tion has been modelled using the Morrison & McCammon 
( 1983| ) cross-sections. The background was taken in a sec¬ 
tor located to NE of G327.1 1.1, with radial distances 
ranging from 10 ' to 20 , (an area reasonably free from point¬ 


like sources), and then corrected for vignetting and extrac¬ 
tion areas. 

We have first extracted a spectrum from a circular re¬ 
gion with 3' radius, centered on the X-ray peak position 
(hereafter referred to as the “plerion” region), with the 
aim of minimizing the contribution of the thermal compo¬ 
nent. For the spectral fitting we have generated an “ad- 
hoc” response file which takes into account the source ex¬ 
tension, that has been estimated in the following way. In 
the 2-10 keV MECS band the apparent size of the source, 
as derived with a gaussian fit, is 3.6' (FWHM), while that 
of the MECS Point spread Function (PSF, obtained by a 
similar fit to the archive observation of Cyg X-l) in the 
same energy range is 2.4': therefore the intrinsic size of the 
source can be estimated as x/3.6 2 — 2.4 2 = 2.7' (FWHM), 
with an uncertainty ~0.3'. 

We have performed different fits (labelled as Fit A, B 
and C) to the spectrum of the plerion region, whose results 
are reported in Table |I]. A pure power-law model using the 
LECS+MECS data (Fit A) provides a good overall fit of 
the spectra in terms of y 2 test, but tends to underestimate 
the flux below 1 keV. The addition of a thermal compo¬ 
nent is needed in order to have the spectra well fitted over 
the whole LECS+MECS bandwidth (Fit B). It is worth 
noting that both the value of the absorption and that of 
the power-law index are significantly different from the 
values obtained with the pure power-law model. The re¬ 
sults of Fit B are confirmed by a pure power-law fit which 
uses only the 2-10 keV data of both instruments (Fit C), 
where the contribution of the thermal component is mi¬ 
nor. Therefore we have decided to take the parameters of 
Fit B as our best guess for the plerion. 

We have then performed a set of fits using a larger ex¬ 
traction region (Fit D and E), that is 6 ' for MECS and 8 ' 
for LECS. The aim of this selection is to include as much 
as possible the contribution of the thermal component. 
The response has been generated consistently, assuming a 
putative size for the shell, estimated using the corrected 
LECS image in the 0.1-1 keV band (Fig. ||-top) and ap¬ 
plying the same procedure as that described for the pleri- 
onic component: The measured, PSF and deconvolved size 
are respectively 21.3', 5.1' and 19.7' FWHM. The results 
of these fits are also reported in Table [3]. Again the two- 
component fit (Fit E) is statistically to be preferred to the 
pure power-law one (Fit D): the y 2 decrease between the 
two is significant at the 99.9% confidence level. The two- 
component fit of this spectrum is shown in Fig. while 
the contour levels for the determination of the tempera¬ 
ture and interstellar absorption are reported in Fig. ^j. The 
contribution of the thermal component to the observed 
spectrum is low. We have estimated that the thermal con¬ 
tribution to the total absorbed 0.5- 10 keV flux is ~ 10%, 
and decreases to 0.4% in the 2-10 keV band. 

When compared with previous X-ray spectral analyses, 
we find a lower value for the temperature and a higher 
value of Nh : Seward et al.| ( 1990| ) in fact derived kT = 
0.8 ± 0.3 keV and TVh = 1-3 ± 0.4 x 10 22 cm -2 , while SWC 
gave kT = 0.371q; 3 o keV an d = 1.8± 0.3 x 10 22 cm -2 . 
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Fig. 2. Folded (left) unfolded (right) spectrum of the of 
G327.1 1.1, with best-fit power-law + thermal models. 
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Fig. 3. Confidence (1-cr, 2 -a and 3-cr level) contours in the 
N]i-T e s parameter plane derived for the total fit (Fit E in 
Table 1) to the LECS+MECS spectrum of G327.1-1.1. 


When the quoted uncertainties are taken into account, 
our results are still marginally consistent with those of 
Seward et al. (1996|) and SWC. The discrepancy should 


be ascribed to differences in the instrumental response of 
LECS (used here) and ROSAT PSPC (used both by SWC 
and by |Seward et al. ( 1996|) for the soft band). 


3.2. The image of the thermal component 

In order to further disentangle the image of the thermal 
component from the total image, in a low-energy image 
we have subtracted the plerionic component, by using the 
following approach. First, we have extracted the MECS 
image in the 2-10 keV energy band, in which the contri¬ 
bution of the thermal component is far below the one of 
the plerion. Then we have convolved that image in order 
to simulate the effect of the LECS PSF at ~ 1 keV (we 
have found appropriate a convolution with a gaussian of 5' 
FWHM). Using the plerion best-fit model (Fit B) reported 
in Table 0, we have scaled its image in order to reproduce 


Fig. 4. 0.1 -1 keV image of the thermal component ob¬ 
tained by subtracting a scaled and smoothed MECS im¬ 
age from the total LECS image in the same band. The 843 
MHz total intensity radio contours of 
( 1996[ ) are overlaid (solid), as well as 
contours at 20%, 40%, 60% and 80% of the peak value. 
The thermal emission extends beyond the plerion in the 
NE and SE direction. 


Whiteoak & Green 
surface brightness 


the LECS image of the plerion only at 0.1-1 keV. Finally, 
we have subtracted it from the original LECS image of 
G327.1-1.1 in the 0.1 1 keV band. The result is shown in 

Fig. 0. 

The image is composed of a large diffuse component 
located in the eastern part of the radio shell and a small 
X-ray blob located near the western rim. The center of the 
large diffuse component is located at R.A.=15 h 55 m 04.1 s , 
Dec.=—55°05 , 08", and it is displaced by ~ 5.5" from the 
radio SNR center and by ~ 5.5" from the plerion center. 
The observed size of this component is ~ 21' FWHM. An 
interesting point is that the image of the thermal compo¬ 
nent does not show a complete shell structure. It could be 
interpreted as partial shell SNR, with a possible extension 
of the X-ray emission in the inner region. The absence of a 
limb can be inferred also from a PSPC image of the whole 
X-ray SNR, as obtained by Seward et al. (199(i|). 


3.3. Timing analysis 

We have searched for pulsation in the events collected in 
the region of the plerion defined above. We have found no 
pulsation at the 99% confidence limit, with an upper limit 
at the same level of confidence of 19.4% in the 0.1-256 Hz 
frequency range. 
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Table 1. Summary of spectral fitting results. 


Name 

Ah 

10 22 cm -2 

7 

0.5-10 keV flux| 

10~ n erg cm -2 s 

kT 

1 keV 

X^/dof 

Fit A - Pie 3' 

Fit B - Pie 3' 

Fit C - Pie 3' 2-10 keV 
Fit D - G327 6' 

Fit E - G327 6' 

‘ >1.4 

O A <3.2 
Z *^>1.7 
o q<3.1 

Z.O>i 7 

U -°>0.6 

2 1 <2,7 

oil <.'Z.'Zia 

Z - 1J ->1.98 
o q't<2.47 
z<z ' >2.08 
o on<3.09 
z - z ^>1.68 

1 Q4 <207 
i - b, ^>1.83 

9 op;< 2 - 44 
z.zu> 2 .07 

1 - 1 ^>0.92 

£O<7100 

oo >0.9 

1 07 <1,57 
x - u ' >0.80 

1 qq<i- 7 o 

i.oy>i is 
/qo< 7 °5 
^>3.0 

U - 1O >0.03 

H 1 o<0-22 
u - i8 >0.14 

139/160 

131/158 

118/141 

235/232 

223/230 


a Unabsorbed non-thermal flux in the 0.5-10 keV band. For Fit A, C, D it refers to the non- 
thermal component; while for Fit B and E we have reported the flux of the thermal component 
only. 


4. Discussion 

4.1. The thermal component 


Table 2. SNR parameters and associated 2<r uncertainties 
derived by the fit to the Sedov and WL model. 


SWC have discussed a scenario in which the shell compo¬ 
nent of G327.1 1.1 is expanding adiabatically in a homo¬ 
geneous medium, and have thoroughly used for it the well 
known Sedov (1959) analytic solution. In this way SWC 
obtained a self-consistent scenario according to which the 
remnant is 1.1 x 10 4 yr old and expands adiabatically in 
an ambient medium with a density of 0.1 cm -3 . 

The two substantially new facts with respect to the 
starting point of their analysis are that the X-ray struc¬ 
ture cannot now be interpreted as a complete shell, and 
that the effective temperature turned out to be lower than 
previously determined. Moreover, with respect to SWC we 
use a slightly higher estimate of the average diameter of 
the SNR, namely 17f5, as derived from the size of the radio 
shell ( |Whiteoak fe Greeii 1996). 

Let us first re-evaluate so me parameters of G327.1- 
1.1 assuming a Bedov ( 1959| ) scenario. In this way we 
derive how the difference from SWC in the starting 
point is going to affect the result. The basic relations 
for the effective temperature (Tea), the age of the SNR 
(tsNR), the emission measure parameter (EM) defined as 
(/ n e rii dV)/{4Ttd 2 ), the radius of the shock (R s h) and the 
total swept mass (M to t) are respectively: 


fcTeff = 1.28* T Bh = 0.14m H R s 2 , (1) 

tsNR = 0.4i? sh /y sh , (2) 

EM = 0.75n^R 3 h /d 2 , (3) 

R sh = 1.15 (E SN fg NR /Po) 1/5 , (4) 

M tot = 4.19 Po R 3 sh , (5) 


where Eg-^ is the supernova energy, no is the atomic num¬ 
ber density in the ambient medium (with standard abun¬ 
dances; po = 1.26wHno), d indicates the SNR distance, 
R s h and T s h are respectively the shock velocity and the 
post-shock temperature: the shock has been modelled as¬ 
suming it to be a strong non-radiative shock with temper¬ 
ature equilibration between ions and electrons. 

Using the information on G327.1-1.1 angular size, the 
above relations translate into the following ones: 

Rsh = 814T^ e y kms -1 , ( 6 ) 


Rsh (kins -1 ) 
tsNR/dio (10 4 yr) 
nod^ 2 (cm -3 ) 
Usid " 572 
M x d - 5/2 (M 0 )| 
d B (kpc) 


Sedov model 


345(304 - 430) 

2.9(2.6 — 3.3) 
0.37(0.10 - 1.35) 
1.4(0.3 - 6.2) 
800(220 - 2910) 

8.8(4.8 - 16.2) 
0.09(0.009 - 1.06) 


WL model 
[C/t = 3.25) 
590(470 - 790) 

1.7(1.2 - 2.1) 
0.11(0.02-0.41) 

3.8(0.7 - 25.1) 
1120(320 - 4000) 

6.0(3.7 — 11.0) 
0.003(0.0001 - 0.14) 


For the Sedov model M x = M to t 


iSNR 

= 12, 200 T k JJ 2 dio£, yr, 

(7) 

n 0 

= 0.511 EM{{ 2 dw / 2 C 3/2 cm" 3 , 

(8) 

Esn 

= 10.70 x 10 51 T keV EMif^/ 2 ^ 3 / 2 erg, 

(9) 

A/tot 

= 1,100 EM Y/d 5 ^ 2 ^ 2 M & , 

(10) 

dE 

= 3.88 £ 5 2 / 5 EM- 1 / 5 T k - e ^ / 5 r 3/5 kpc, 

(11) 


which directly link physical and observable quantities. 
In these equations T ke v is the effective temperature, in 
keV; EM 14 is the emission measure parameter, in units 
of 10 14 cm -5 ; dio is the G327.1-1.1 distance, in units of 
10 kpc; and £51 is the supernova energy, when assumed, 
in units of 10 51 erg. The quantity £ indicates the SNR an¬ 
gular size, in units of the G327.1 1.1 one. The physical 
quantities as derived from our spectral fit in Sedov ap¬ 
proximation are listed in the first column of Table |2j. It is 
apparent that, while the swept mass is much larger than 
any reasonable stellar mass, the SNR shock is slow, so 
that the resulting supernova energy is rather normal for 
any reasonable distance value. 

In other terms, assuming the standard value of 10 51 erg 
for the supernova energy, we derive a distance between 4.8 
and 16.2 kpc, which is in agreement with published dis¬ 
tance estimates, ranging from 6.5 kpc (Seward et al. 1996) 
to 9 kpc (SWC), as well with the distance estimated from 
the T.-D empirical relation ( |Case fe Bhattacharya 1998). 
The average 1 GHz surface brightness of the G327.1 1.1 
shell component is 2.4 x 10 -21 W m -2 Hz -1 sr _1 (5.6 Jy at 
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T. ff (keV) 


Fig. 5. Efficiency of radiative cooling (e c , scaled with 
d 1 o ), evaluated using Sedov solutions. The solid straight 
lines show levels of the radiative efficiency in the Tk e v 
EM plane. The best fit solution (T e g- = 0.18 keV; EM = 
0.52 x 10 14 cm -5 ) is indicated by a cross. The confidence 
contours are at 1-cr, 2-cr and 3 -a level (computed from Fit 
E in Table g). 


843 GHz, with a spectral index ~ —0.4; SWC): this trans¬ 
lates into a shell diameter of 45 ± 20 pc, and then into a 
distance estimate of 8.8 ± 4.0 kpc. 

To summarize the results of our Sedov fit compared to 
those of SWC (using the same assumed distance of 9 kpc), 
we find a SNR ~2.4 times older and with ~4.7 times more 
energy, expanding in a medium ~3.9 times denser. The 
swept mass is 12.5 times larger than that derived by SWC. 

The last row in Table g gives the efficiency in radiative 
cooling, defined as the ratio between the SNR age and the 
present cooling time right behind the shock, which can be 
expressed in terms of the observed quantities as: 

e c = 2.8 x 10 _3 (T k ((y 2 + 0.33T k ^) v /EM 14 cW£, ( 12 ) 


where for the cooling coefficient we have used the formula 
A (T 6 ) = 1.0 x 10” 22 T 6 "°- 7 + 2.3 x 10 


24 T °- 5 erg cm 3 s 4 , 


with Tfj in units of 10 6 K ( McCray 1987 ). A value for e c 
much smaller than unity is required in order to test the 
adiabaticity condition, as required by the Sedov model. 
However Fig. g shows that, since e c is strongly decreasing 
with T, fits giving low values for T e g also imply relevant 
radiative losses. 

The possible extension of the thermal X-ray compo¬ 
nent to the inner shell regions, together with the non- 
negligible radiative cooling, leads us to consider scenarios 
alternative to the Sedov one for the evolution of G327.1- 
1 . 1 . 


Let us first examine the WL model, describing (by a 
class of self-similar solutions) a wide range of conditions 
under which a SNR may expand in a cloudy interstellar 
medium. In this scenario, a delayed evaporation of the 
clouds causes the central density to increase with respect 
to the pure Sedov case. This class of solutions depends 


on two parameters: the mass fraction in clouds ( C ) and 
the cloud evaporation time scaled with the SNR age (r). 
The solutions converge to the Sedov case when the mass 
deposition is small, either because little mass is contained 
in the clouds (C <C 1) or because the clouds have anyway 
released little mass during the SNR lifetime (C <C r). 
Here we assume that the mass in clouds is large (C> 1), 
but the evaporation time is long (r /§> 1 ): in this case 
the solutions form a one-parameter class, depending on 
just C/t (these solutions are approximately valid also for 
C,tZ 1 ). 

We have qualitatively compared the radial profile of 
the thermal component (using the center coordinates re¬ 
ported in Sect. 3.2) with the WL predicted profiles for 
various value of the C/t parameter. We found a reason¬ 
able matching, considering the uncertainties involved, for 
C/t in the range between 2.5 and 4.0. These uncertain¬ 
ties originate not only from the faintness of the source, but 
also from deviations from the assumed spherical symmetry 
as well as from a non perfect subtraction of the plerionic 
component: all of these effects, with the present data, are 
even difficult to quantify. Therefore, under the assumption 
that WL models are appropriate to model the G327.1 1.1 
thermal component, we can estimate C/t = 3.25 ± 0.75. 
By using the results tabulated by WL we derive which 
correction factors must be applied to the Sedov estimates. 

For instance 

R sh = R sh /R sh , s = (0.81 ± 0.04) (13) 


shows how much the shock radius shrinks with respect to 
the Sedov case (for fixed Esn, no and fsNR) : with reference 
to the quantity K, defined in Eq. 7 of WL, the relation 
i? 3 h = K/K$ holds. Moreover the equations 

EM = EM/EMs = (10 ± 8 ); (14) 

T eff = T eff /T eff , s = (0.33 ± 0 . 12 ); (15) 

M x = M x /M totis = (4.5 ± 1 . 6 ) (16) 


refer to comparisons with the Sedov solution keeping no 
and i? s h constant. The above scaling factors are evaluated 
by interpolating data from Table 6 of WL. In the above 
expressions we have used the index g to indicate the values 
in the Sedov solution. 

Therefore the corrections required to derive, for the 
WL model, the relations equivalent to Eqs. 6-11 for the 


associated Sedov case are: 

H sh = ?eff 1/2 ^h,s = (1-7 ± 0.3) Fsh.s, (17) 

^SNR = T e 1 f £ 2 tsNR,S = (0.6 ± 0.1) tsNR.S, (18) 

-— 1/2 

no = EM n 0 ,s = (0.3 ± 0.1) n 0j s, (19) 

£ sn = R~ E EM _ 1/2 T~g Esn,s = (2.7 ± 0.7) £ SN , S , (20) 
-- — 1/2 

M x = M x EM Aftot,s = (L40 ± 0.04) M tot , s , (21) 
d E = i? 2 h EM 1/5 f e 2 ff /5 d E ,s = (0.68 ± 0.07) d E ,s : (22) 


It should be noted that here by “associated” WL and 
Sedov solutions we indicate those showing the same ob¬ 
served properties (i.e. size, T e ff and EM). 
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The values of the parameters derived with the WL 
model are also given in the second column of Table |j. 
Compared to the pure Sedov case, now the ambient (inter¬ 
cloud) density results to be lower, while the shock velocity 
is higher (which also makes the SNR younger). The dis¬ 
tance estimate based on the 10 51 erg energy requirement 
now gives ~ 6.0 kpc, still compatible with the lowest dis¬ 
tance estimates for this source ( {Seward et al. 1996| ). The 
cooling efficiency e c (in the intercloud medium) can be 
evaluated from Eq. 12, by substituting Tk e v with Tkev/Teff 
and EM 44 with EM 44 /EM: it results much smaller (~ 
3 x 10 -2 ) than in the Sedov case. 


On the basis of our analysis, a WL model fits the aver¬ 
age radial profile of the Thermal remnant better than the 
Sedov model. Moreover the WL model gives reasonable 
values of the SNR parameters, even though: 1. it leads to 
a rather small SNR age, with problems for the offset of 
the associated neutron star (see below); 2 . for a standard 
SNR energy, it leads to a distance somehow smaller than 
what estimated by other methods. 

A further scenario which may account for thermal X- 
ray from inner regions is that of a SNR in radiative phase 
(with in addition possible effects, like thermal conduction 
and flux saturation). A common property of non-radiative 
spherical shocks expanding in a homogeneous medium is 
the presence in the inner region of a hot and thin medium, 
which therefore behaves as a very low efficiency emitter, 
and forms the typical hollow sphere emission pattern of 
shell-type SNRs. In a cloud evaporation model (as WL, see 
above), this region is filled with material originally stored 
in small clouds. But since inside the shell the pressure 
equilibrium is roughly maintained, it would be sufficient to 
cool down this gas (even without adding further material) 
in order to make it denser. 


Unfortunately no self-similar solution is known for ra¬ 
diative shell-type SNRs, and therefore a quantitative de¬ 
scription of this class of objects requires to develop nu¬ 
merical models, which are beyond the scope of this work. 
Models of this kind have been already developed, for in¬ 
stance for W44, by various authors (Harms et al. 1997; 
Cox et al. 1999; |5helton et al. 1999| ). Here let us simply 
compare the basic parameters of G327.1 1.1 and W44 to 
infer their relative evolutive conditions, under the assump¬ 
tion that both have reached the radiative phase. 

Table ^ lists, for both remnants, some basic parame¬ 
ters, either observed or derived: in the latter case a Sedov 
model has been used for the derivation. Thus if these SNRs 
are no longer in Sedov phase, the listed values are inaccu¬ 
rate (even though, simply for dimensional arguments, not 
too far from the correct values). 

For instance the absolute ages, as reported in Table ||, 
may be understimated since in radiative phase the SNRs 
must have decelerated with respect to the Sedov phase. 
We know this is the case for W44, where the Sedov es¬ 
timate is about half of what derived from more detailed 
NEI models, as well as half of the pulsar spin-down age 
(Harms et al. 1997). 


Table 3. Comparison between G327.1-1.1 and W44 (from 


model); for W44 we used Eqs. 

6 through 12 with £ = 


G327.1 1.1 

W44 

Ash 

25.5 d 10 pc 

11.2 d 2 .5 pc 

kT 

0.18 keV 

0.88 keV 

{n 2 )V 

7.0 x 10 59 di 0 cm -3 

1.3 x 10 58 d|. 5 cm -3 

ilGHz 

6.2 x 10 32 di 0 ergs -1 

1.7 x 10 33 d!. 5 ergs -1 

age 

2.9 x 10 4 di 0 yr 

5.7 x 10 3 d 2 . 5 yr 

no 

0.37 d~g /2 cm -3 

0 . 18 d^ 5^ 2 cm -3 

e c 

0.09d^ 2 

7.5 x 10 -4 d 1/2 


If also G327.1 1.1 is in radiative phase, its true age 
may be substantially larger than 29,000 yr (as from 
Table ||). At any rate, if we exclude the WL scenario, 
this SNR must be much older than originally estimated 
by SWC (11,000 yr). An implication would be that the 
lower limit on the velocity of the compact object (> 600 
km s -1 , as from SWC, in order to account for the displace¬ 
ment of the compact source from the SNR center) can be 
lowered by a factor 3 or even more: therefore in this case 
there is no need for an anomalous pulsar velocity. 

Although some values listed in Table |] may be inac¬ 
curate in an absolute sense, we are confident that they 
may still be used for a comparative analysis of the two 
SNRs. For instance G327.1 1.1 looks older than W44; it 
also expands with a lower velocity, and in a slightly thin¬ 
ner medium. Moreover, as shown in the last line of the 
table, the cooling efficiency in G327.1-1.1 is much higher 
than in W44, and this fact strengthens the idea that, if a 
radiative scenario results to be appropriate to the case of 
W44, it should be even more reasonable for G327.1 1.1. 

Although expanding in a higher ambient density, 
G327.1-1.1 radio shell is about 3 times fainter than W44: 
also this difference could be a consequence of its age. The 
fact that high resolution maps of the radio continuum from 
W44 present a strongly hlamented pattern has been taken 
( pox et ah 199E) as evidence for the remnant being in ra¬ 
diative phase. If this is the case, higher resolution radio 
maps of G327.1 1.1, able to test the presence of filamen- 
tation, could be valuable to shed light on the evolutive 
phase of this remnant. 


4.2. The plerionic component 

Like W44, G327.1-1.1 shows both a shell and a plerionic 
component in radio, while a thermal filled-centre and a 
non-thermal component in X rays: the main difference is 
that in G327.1 1.1 the plerionic component is more promi¬ 
nent than in W44. Another similarity between W44 and 
G327.1-1.1 is the displacement of the plerionic component 
from the shell centroid, as well as its bow shock shape. 

An analysis of the integrated spectrum of the plerionic 
component can be carried on in a very similar way as 
already done by SWC. Our estimate of the 0.5-10 keV 
luminosity of the plerion (~ 1.6 x 10 35 df o ergs -1 ) is about 
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0.5 times that estimated by SWC; moreover the age we 
find is at least a factor 2.4 larger, and all this affects the 
derivation of pulsar and nebular parameters. 


By applying an empirical relation (Seward & Wang 


1988), SWC used the X-ray nebular luminosity (Lx) to 


infer the pulsar spin down power (E), which for a magnetic 
dipole braking is proportional to P/P 3 . The age estimate 
(^SNR ~ P/2P) is the other relation which allows one 
separating the pulsar timing parameters P and P, and to 
estimate also the pulsar surface field: SWC report P = 
62 ms and B 0 = 2.3 x 10 12 G. Since by this procedure P oc 
fsNR 2 ^ -1 / 2 ; while Bo oc f S jJ R P 1 / 2 , our revised values 
lead to a somehow shorter period (< 50 ms) and to a lower 
surface field (< 1.2 x 10 12 G). The presence of a faster, 
and lower-B pulsar compared to that in W44 is in fact 
required by the presence of a reasonably bright, although 
aged, plerion. 

SWC also used a spectral break in the nebular spec¬ 
trum (at about 3.5 x 10 4 GHz) to infer a nebular field 
B ~ 0.7 x 10~ 4 G. Since the field derived in this way is 


oc t. 


-2/3 
SNR i 


we infer instead a nebular field ^ 0.4 x 10 4 G. 


For a radius of ~ 1.5arcmin for the (radio) plerion, the 
total magnetic energy is ^ 6 x 10 4 ' erg. The plerion is 
then slightly underpressured (Pb < 0.2x 10 -10 dyncm^ 2 ), 
compared to the thermal remnant, if in Sedov or WL 
regime; while a better pressure balance may be attained if 
the shell is in radiative expansion. Anyway the uncertain¬ 
ties involved in these estimates are too large to take this 
as a strong argument in favour of a radiative shell. 


5. Summary and conclusions 

We have used a BeppoSAX X-ray observation of G327.1- 
1.1 to study the origin of its X-ray emission, with partic¬ 
ular emphasis on its thermal component, which has not 
been studied in detail in the literature, because of its su¬ 
perposition with the bright non-thermal plerion. We fitted 
the plerion and the total SNR spectrum with a combina¬ 
tion of thermal and non-thermal models, finding an X- 
ray temperature kT ~ 0.18 keV for the hot plasma and 
a power-law photon index ~ 2.2 for the plerion. Using 
the flux derived for the non-thermal plerion, we have sub¬ 
tracted the non-thermal component from the total 0.1-1 
keV SNR image, thus obtaining for the first time an image 
of the thermal hot plasma. 

The emission from the plasma is not confined to a 
bright shell, and shows possible emission for inner regions. 
We have then examined for this SNR also scenarios differ¬ 
ent from the Sedov one. The model of SNR expansion in 


a medium w ith evaporative clouds developed by White fe 
Long (1991) may account for the observations, although it 
may imply an age and distance for this object which are 
somehow smaller than expected. We have shown that a 
radiative expansion with strong central cooling may also 
describe the data. We have derived and compared the 
remnant characteristic parameters for the various scenar¬ 
ios, either Sedov, or White fe Long| (1991), or radiative 
model, showing that the last model implies a very long 


age (ri 3 x 10 4 yr): in this case the observed offset of the 
associated pulsar may be accounted for without requiring 
a too high spatial velocity for it. 

A comparison with W44, a remnant for which there 
is also evidence of radiative expansion, indicates that 
G327.1-1.1 may be substantially more evolved than W44. 
Future high sensitivity X-ray and radio observations are 
strongly encouraged, because they may dissipate any 
doubt about the evolutionary stage of this object and in 
particular of its shell. Finally, we have also revised the 
estimate of the pulsar and nebular magnetic field. 
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